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PPJiFACE 



GacPl i^H t'ccomes c?,>f fleering problem. The ^ i - 't •jrdi.n ii i ly srist;s 
in the ge/ir crain connect the output fjotoc " • Isad*, The fer- 

foriiinuct^ .changes resulting from the addition f f ’ nrl^^.rsh depend upon 
whether, i he najor portion ot the output inertia ii on the motor or tlie 
load zLle. of the backlast^. Vo matter where the c«<:*/Ginc chan ism error is 
T!if*?^S‘ire 1, h'>rklash contrihutei5 eirher to the i>rrrdy r;tate error^, or in- 
sr<*blltf: the system. 

BacLlanh effects are investigated by analog coi;»p.ur-. r. techniques in this 
tbasu s o 

We wish to acknowledge the assistance, the corim.cnt s ^rd suggestions we 
have received from ProfesJot Oto George J« Tlmiti that make us 
possible to complete rhis thesis^ 



il 



1 







TABLE OF CONTENTS 



Section pggg 

1. Background o£ the problem. 1 

2. Review of second order system with backlash 14 

and load inertia. 

3. Physical analysis of third order system with 27 

backlash and load inertia. 

4. Analog computer technique. 65 

5. The results of investigation. 75 

6. Discussion and conclusions. 78 



ill 



List of Illustrat ions 



Page 

1. Fig. 1 Block diagram of a simple feedback control system 3 

2. Fig. 2 Relationship between angular position of two gears 3 

3. Fig. 3 A second order servo with backlash 5 

4. Fig. 4 Phase plane trajectory of a second order servo 6 

5. Fig. 5 A sketch of displacement vs. time 7 

6. Fig. 6 A sketch of velocity vs. displacement 10 

7. Fig. 7 A sketch of acceleration vs. displacement 11 

8. Fig. 8 A sketch of tofque vs. displacement 12 

9. Fig. 9 A sketch of displacement vs. displacement 13 

10. Table 1 18 

11. Table 2 19 

12. Table 3 20 

13. Table 4 23 

14. Fig. 10 A graph of displacement vs. time 25 

15. Fig. 11 A phase plane trajectory of velocity vs. displace- 26 

ment. 

16. Fig. 12 A block diagram of a third order servo 29 

17. Fig. 13 A phase plane trajectory vs. displacement 33 

18. Fig. 14 A plot of acceleration vs. displacement 35 

19. Fig. 15 A plot of motor velocity vs. displacement 41 

20. Fig. 16 A plot of load displacement vs. motor displacement 42 

21. Fig. 17 A plot of motor torque vs. motor displacement 43 

22. Fig. 18 The root locus of the system 44 

23. Table 5 48 

24. Table 6 52 

iv 



List of illustrations (Cont'd) 



Page 

25. Table 7 53 

26. Fig. 19 A complete phase plane trajectory with backlash 54 

.10 rad, 

27. Fig, 20 A complete phase plane trajectory with backlash 55 

.20 rad. 

28. Fig. 21 A complete phase plane trajectory with backlash 56 

.30 rad. 

29. Fig. 22 A complete phase plane trajectory with backlash 57 

.40 rad, 

30. Fig, 23 A complete phase plane trajectory with backlash 58 

,50 rad, 

31. Fig, 24 A complete phase plane trajectory with backlash 59 

.60 rad. 

32. Fig, 25 A complete phase plane trajectory with backlash 60 

.70 rad. 

33. Fig. 26 A complete phase plane trajectory with backlash 61 

.80 rad. 

34. Fig. 27 A complete phase plane trajectory with backlash 

.90 rad. 

35. Fig. 28 A complete phase plane trajectory with backlash 63 

1.0 rad. 

36. Fig. 29 Trasient response of motor and load after separation. 64 



37. Fig. 30 A circuit diagram of computer set up for the 67 

combined system. 

38, Fig. 31 A circuit diagram of computer set up for the 72 

load and motor after separation 

39. Table 8 76 

40, Fig, 32 The plot of amplitude of limit cycle vs. 77 

backlash with several value of gain 



V 



lo Background of problem« 

In the study of backlash we have to consider two cases? 

a) Simple backlash theory, in which there is no inertia on 
load side« 

b) Complicated backlash i.e« inertia and friction on load 

sideo 

ao Simple backlash theoryo 

A simple feedback control system is assumed as shown in Fig. 1. 
Backlash occurs in gear train connecting the output motor to the load. 

The summer is assumed ideaT here in its function of sending to the con- 
troller the difference between input signal and output of box containing 
gears with backlash. The input to the system is a shaft rotation, de- 
signated in Fig. 1. The output shaft is a rotation and designated 

with -0 . The output of the gear box, whose output is O , is assumed 

t» V* 

ideal in all respect except for the presence of backlash. Because there 

is no inertia on load side the backlash is such a character as shown in 

Fig. 2. This figure shows that as the output shaft moves with finite 

velocity from a position 0 is zero to a position 0 has some posi- 

tive value, say greater than 1/2 the backlash angle, the output of 

the gear box 0® moves from zero to a certain positive value; then the 

direction of rotation reverses so that 0^ decreases through zero to 

some negative value than 1/2 the backlash angle. The value 

of 0® goes from positive through zero to negative value and then goes 

back to zero again. is the total backlash angle. There is a dead 

space and this means that the adjacent gear teeth must move through a 

finite distance before making contact with each other. As long as the 

operation is unidirection there is no trouble, but soon when it changes 

its direction, then backlash must be taken up each time. In manual 

1 
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operation it is sometimes noisy and can cause instability in the systemo 
bo Qualitative effect of inertia on loado 

Assumed vje have a unity feedback system as shovm in Figo lo 
The discussion of the physical nature of the system is as follows! 

In normal operation the motor may be driving in one direction and 
backlash is then taken up^ so that all gears are in contact. If the load 
inertia is small and there is some load friction, the load shaft will stop 
almost immediately with the reversal of the motor rotation. The load shaft 
will remain stationary until the motor has moved through an angle suffi- 
cient to take up all backlash. On the other hand, if the load inertia is 
big enough, and the load friction is small, when the motor reverses the 
load will not reverse, but may continue to rotate in its original direc- 
tion with a constant velocity equal to its velocity at the instant the 
gear tooth contact was broken. If both the load inertia and friction are 
appreciable, the load rotates in its original direction but with a con- 
stant deceleration. Under these conditions the motor shaft and load shaft 
both help to take up the backlash. 

c. Second order servo with backlash. 

The schematical diagram of a secondorder servo is shown in Fig. 

3. In the normal arrangement the load is geared to the motor through 
gear train and the load shaft is geared to a measuring device 

through gear train G 2 . There may be backlash in either of these gear 
trains or in both of them. For our purpose we consider the backlash in 

Gr. only. Consider the sketch of Fig. 4. Assuming a step displacement of 
2 

magnitude OA and also assuming that the backlash happens to be all taken 
up in the proper direction, the trajectory starts at A. The trajectory 
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Fig. 2. Relationship between angular position of two gears between 
which backlash exists. 
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follows a linear path. At point B the velocities of the motor and 
load become different, mechanically contact is lost and separation 
occurs, providing the friction of load has a finite value and the coloumb 
friction of load can be neglected. The point of separation is determined 
by the load trajectory, when this trajectory has a common slope with the 
trajectory of the system. Combining these points together it forms a 
straight line and is called the "Separation dividing line." After 
point B the system operates in the backlash region. If the viscous 
friction of the load is negligible the trajectory of the load is a 
straight line with constant- velocity. But when the friction is suffici- 
ent enough the load is drifting separately with a constant deceleration. 
The remainder of the system continues to drive in the reverse direction 
until the backlash has taken up and the system output completely cor- 
rected as to displacement. The separation dividing line determines the 
initial conditions for load velocity and displacement, which are the same 
value as the initial conditions for the motor velocity and displacement. 
The sketch of Fig. 5 shows the displacement of load as a function 
of time and also the motor displacement as a function of time. From 
these two curves the time at which the backlash is taken up can be deter- 
mined. Knowing this time the velocity and displacement of load at time 
of recombination can be determined. This is point C in the sketch of 
Fig. 4. At the same time point D is also located. This point Indicates 
the velocity and displacement of the system without load. E 4 (ch pair of 
corresponding points for load and system without load, velocity and dis- 
placement at recombination, in turn determine the velocity and displace- 
ment for the recombined system. This is obtained by applying the law of 
conservation of momentum. The point of recombination for the system is 
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a 3* A second order servo with backlash 
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Load- d i sp lacement 




Fig. 5. A sketch of curves of displacement vs time. 
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point E in Fig* 4* Trajectories for the recombined system start from 

E and the system retains its original characteristics until the next 

separation dividing line is reached; separation and recombination will 
occur again as in the previous half cycle. If the final state of the 
system spiral into the origin the system is said to be stable. But when 
the final state of the system oscillates with a constant amplitude, the 
system is said to have a limit cycle. 

d. Higher order system. 

For higher order, because of the complex mathematics we have 
to break it up into severaT. phase plane trajectories in order to analyze 
the system. For example: 

1. The plot of velocity vs. displacement (Fig. 6) 

2. The plot of acceleration vs. displacement (Fig. 7) 

3. The plot of torque vs. displacement (Fig. 8) 

The conditions at the point of separation is the same as the second 
order system. At the point of recombination, the law of conversation 
of mementum is still applicable, but we also have to consider the torque 
equation of the motor, In order to determine the acceleration. Knowing 
the acceleration and velocity the point of recombination is located. The 
trajectory of the combined system start again until the separation divid- 
ing line is reached. The picture repeats again as in the previous half 
cycle. 

Note here, that for the determination of the load displacement and 
velocity it is not necessary to go back to the time function. In order 
to know the value of load displacement and motor displacement at which the 
backlash is taken up, just plot the load displacement versus motor displace 

ment as shown in Fig. 19. Using this value we go back to Fig. 6, where 
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the velocity of motor is plotted against the motor displacement* For 
the calculation of velocity at recombination, the law of conservation 
of momentum is applicable as the case in the second order system* At 
the recombination point the torque equation of the motor is true, since 
the load has no torque at all and also providing the current of the 
motor is a continuous function of time during the separation. From 
Fig. 8, we can determine the value of the torque at the time of re- 
combination. Knowing the torque and the velocity, we can evaluate the 
value of acceleration at this point. 
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velocity 
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Figo 7o A sketch of acceleration vs, displacement of 
the combined system. 
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torque 




> 

time 



Flgo 8o A sketch of torque vs« displacement of the motor. 
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load displacement 




Flgo 9o A sketch of load displacement vs<> motor 
displacement 
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2o Review of Second Order System with Backlash and Load Inertia* 
For our particular problem assume the gear ratio is equal to 
the equation for the combined system becomes: 



(■'m * \ > \ * * 'l> -®C + = >">R 



For the system with load removed: 



M M 






K6 



M 



= K6, 



For the load drifting separately: 









0 



Given: + J, = 1 + f, = 0.40 K = 1.0 

ML ML 



J„ = 0.50 = 0.30 •dj, = 1.0 

M MR 

J = 0.50 = 0,10 A = 0.30 

L L 



Equation (1) becomes; 

le^ + 0.400^, + 1 



To obtain isoclines for the combined system; 



e. 




+ 0.40 




one. 



( 1 ) 



( 2 ) 



(3) 



rad. 



(4) 



(5) 
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Let N 



1 



N, + 0.40 = 



1 - e. 



(5a) 



Computed values of 6 and 0 for various values of N. are listed 

L/C./ J, 

in Table one. 



The system as open loop with load separated: 



0.500^ + 0.300^ = 1 - e, 

M M L 



( 6 ) 



At the instant of separation; 



e, 



M + 0.30 



1 - 0 



Lo 



-0 



M 






Mo 



(7) 



At the instant of separations 



e ==0 = e 

Co Mo Lo 



^Co = % 



= 0 



Lo 



Hence the equation (7) becomes; 



N., + 0.60 = ' ^Mo^ 



0 



Mo 



(7a) 



where N. 



0, 



il 

M 



15 



2 



^Mo 

1 - -H 0,60 



( 8 ) 



Upon separation of the load, equation (3) becomes: 



0.506- + 0,100- = 0 

L 



(9) 



N, 



L 

T 



= - 0.20 



( 10 ) 



Letting N 2 = == -0.20, equation (8) becomes; 



0 , 



MO 



== 5.0 



1 - 0 , 



Mo 



hence 0„ = ( 1 - 0.20) 6 

Mo Mo 



( 11 ) 



arc, tan. Mo 



= 78.7' 



Mo 



/ = 101 . 3 ‘ 



Solve for equation )9) as a function of time we obtain: 



0.50 ( s 0, 



s0, - 0, ) + 0.10 ( s0, - 0, ) = 0 (12) 

Lo Lo L Lo 



s ( .5 s + .1 ) 0j^ s ( .5 s + .l)»Lo ^Lo 



where 0_ = 0., = 0„ 

Lo Mo Co 

Equation (13) becomes: 

1 



0. 



Mo 



0 , 



s ( s - .2) 



Mo 



(14) 
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Substitute the value of equation (11) into equation (14) it becomes: 



hence : 



Values 



e. 



( 1 - .2 ) 



©M + 

Mo 



€ 



s ( s - .2) 



Mo 



( 15 ) 



0^ = (1 - .20) L + 5.0 (1 - 

L Mo Mo 



(16) 



of versus time (t) are tabulated in Table 1. 
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TABLE 1. 



t 


--2t. 
(1 - e ) 


5(1 - e'*^'') 


-.2 




.2 


.0392 


. 196 


-.2 


.004 


.4 


.077 


.385 


-.2 


.185 


.6 


.113 


.565 


-.2 


.365 


1.0 


.181 


.905 


-.2 


.705 


1.4 


.244 


’l.24 


-.2 


1.04 


1.8 


.302 


1.52 


-.2 


1.32 


2.2 


.356 


1.78 


-.2 


1.58 


2.6 


.405 


2.205 


-.2 


1.825 


3.0 


.450 


2.25 


-.2 


2.05 


3.4 


.493 


2.465 


-.2 


2.265 


3.8 


.532 


2.66 


-.2 


2.46 


4.0 


.550 


2.75 


-.2 


2.55 



Fig. 10 is a graph of equations (16) and (17) . By appropriate ordinate 
scaling, recombination times were computed as listed in table number four; for 
corresponding values of V.alues of for dividing lines 

are listed in Table four. 

Substitute the value of equation (15) into equation (6) and solve for ©^^ 
as a function of time: 

©,, = - 16 ©,, t + 111,455 6^ + 13.345 ©' - 125 0 e'’^*'+ 1.00 

M Mo Mo Mo Mo 

(17) 

Values of 6 and 0., versus time (t) are tabulated in Table three. 

M Mo 
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TABLE 2 



t 


-. 2 t 

e 


-.65 

e 


-125 

e 


13.345 


- 16 t 


111.455 




.2 


.9608 


.887 


-120 


11.82 


- 3.2 


111.455 


.075 


.A 


.923 


.787 


- 115.15 


10.5 


- 6.4 


111.455 


.105 


.6 


.887 


.698 


- 110.9 


9.305 


- 9.6 


111.455 


.260 


.8 


.8525 


.620 


- 106.5 


8.28 


- 12.8 


111.455 


.430 


1.0 


.819 


.550 


- 102.1 


7.35 


-16 


111.455 


.70 


1.2 


.787 


.487 


- 98.25 


6.50 


- 19.2 


111.455 


.505 


l.A 


.756 


.432 


- 94.5 


5.76 


- 22.4 


111.455 


.315 


1.6 


.726 


.384 


- 90.85 


5.12 


- 25.6 


111.455 


.125 


1.8 


.698 


.340 


- 87.25 


4.545 


- 28.8 


111.455 


-.05 


2.0 


.670 


.302 


- 83.8 


4.035 


- 32.0 


111.455 


-.310 


2.2 


.644 


.268 


- 80.45 


3.58 


- 35.2 


111.455 


-.615 


2.4 


.620 


.237 


- 77.5 


3.16 


- 38.4 


111.455 


- 1.285 


2.6 


.595 


.210 


- 74.45 


2.805 


- 41.6 


111.455 


- 1.790 


2.8 


.572 


.187 


- 71.5 


2.50 


- 44.8 


111.455 


- 2.345 


3.0 


. 550 ' 


.166 


- 68.8 


2.22 


- 48.0 


111.455 


- 3.125 


3.2 


.528 


.147 


- 66.0 


1.961 


- 51.2 


111.455 


- 3.784 


3.4 


.507 


.130 


- 63.4 


1.735 


- 54.4 


111.455 


- 4.61 


3.6 


.487 


.116 


- 60.95 


1.55 


- 57.6 


111.455 


- 5.545 


3.8 


.468 


.102 


- 58.5 


1.361 


- 60.8 


111.455 


- 6.484 


4.0 


.450 


.091 


- 56.2 


1.213 


- 64.0 


111.455 


- 7.532 
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TABLE 3 



t : 


_ -.25 

25 e 


o • 6 1 

- 8 e 


- 16 




.2 


24 


- 7.1 


-16 


.90 


.4 


23.03 


- 6.30 


-16 


.73 


0 6 


22 . 15 


- 5.575 


-16 


.573 


.8 


21.3 


- 4.96 


-16 


.34 


1.0 


20,42 


- 4.40 


-16 


.02 


1.2 


19.65 


- 3.895 


-16 


- .245 


1.4 


18.90 


- 3.455 


-16 


- .445 


1.6 


18. 15 


- 3.065 


-16 


- .915 


1.8 


17.41 


- 2.72 


-16 


- 1.31 


2.0 


16.73 


- 2.42 


-16 


- 1.69 


2.2 


16.07 


- 2.14 


-16 


- 2.07 


2.4 


15.49 


- 1.895 


-16 


- 2.405 


2.6 


14.86 


- 1.68 


-16 


- 2.82 


2.8 


14.30 


- 1.495 


-16 


- 3.195 


3.0 


13.73 


- 1.33 


-16 


- 3.60 


3.2 


13.19 


- 1.18 


-16 


- 3.99 


5.4 


12.66 


- 1.04 


-16 


- 4.38 


3.6 


12.16 


- .929 


-16 


- 4.769 


3.8 


11.69 


- .816 


-16 


- 5.126 


4.0 


11.22 


- ,728 


-16 


- 5.508 
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TABLE 3 



^Mo 


t 


“o2t 

e 


- . 6 t 
e 


-. 25 e "‘ 


2 t 

13.345 


-. 6 t ^ 

e - 16 t 


111.455 




1.0 


1.1 


.803 


.517 


- 100.4 


6.9 


- 17.6 


111.455 


.355 


.8 


1.175 


.791 


.494 


- 99.0 


6.6 


- 18.8 


111.455 


.204 


.6 


1.275 


.775 


. 466 


- 96.9 


6.225 


- 20.4 


111.455 


.228 


.4 


1 . A 25 


.752 


.426 


- 94.0 


5.7 


- 22.8 


111.455 


.142 


.2 


1.875 


.688 


.3245 


- 86.0 


4.335 


- 30.0 


111.455 


.044 


.1 


2,425 


.616 


.234 


- 77.0 


3.125 


- 38.8 


111.455 


-.122 


.08 


2,64 


.590 


.206 


- 73.8 


2.75 


- 42.25 


111.455 


-.1475 


.06 


3.2 


.528 


.147 


- 66.0 


1.96 


- 51.2 


111.455 


-.227 


.04 


3.525 


.494 


.1205 


- 61.6 


1.61 


- 56.5 


111.455 


-.201 


.03 


3.98 


.452 


.092 


- 56.5 


1.23 


- 63.6 


111.455 


-.2225 
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TABLE 3 continued 



e „ 25 - 8.0 e - 16 

Ho M 



0 


20.08 


- 4.135 


-16 


-.055 


8 


19.8 


- 3.94 


-16 


-.112 


1 6 


19.35 


- 3.73 


-16 


-.222 


4 


18.18 


- 3.41 


-16 


-.244 


,2 


17.2 


- 2.60 


-16 


-.28 


,1 


15.4 


- 1.87 


-16 


-.247 


08 


14.75 


1 65 


-16 


-.232 


06 


13.2 


- 1.175 


-16 


-.238 


04 


12.32 


-.965 


-16 


-.151 


03 


11.3 


-.735 


-16 


-.113 
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TABLE 3 continued 



°Mo 


t 


KM 




1.0 


1.1 


.803 


.785 


.8 


1.175 


.632 


.675 


.6 


1.275 


.465 


.555 


.4 


1.425 


.30 


.4155 


.2 


1.875 


.1375 


.272 


.1 


2.425 


, .0616 


.172 


.08 


2.64 


.047 


.148 


.06 


3.2 


.0316 


.1292 


.04 


3.525 


.01975 


.0933 


.03 


3.98 


.01352 


.0761 


To satisfy 


the principle of 


conservation of momentum, 


applied the next 


equation; 




■’l.> - Vm + Vl 


(18) 



From equation (18) 6^^ .and 0^ are known, hence we can solve this equa- 
tion. The values of 6^ upon recombination corresponding to values of 
^Mo separation are listed in Table l\ 

TABLE 4 






e, 



Mo 



c 



1.0 

.8 

.6 

.4 

.2 

.1 

.08 

.06 

.04 

.03 



.374 

.260 

.1216 

.029 

.07125 

-.0927 

-.0925 

-.1032 

-.0656 

-.0497 
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Fig. tl is the phase plane presentation of 0,, 6 , and 0 

*L> Li 

as a conbincu system and when operating in the backlp^sh region, A step 
input of 0o80 is seen to separate the load from the system at A, recombina- 
tion takes place when the load is at B, and the balancing of momentum 

« 

between load and the system without load, causes ^ , 0^ of the re- 

combined system to originate from point C. Reseparation occurs at point 
D and the trajectory is seen to spiral into a limit cycle defined by 
points E, P, G and H. The magnitude of the limit cycle is 0.550. 
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Fig, Ho A phase plane trajectory of 0 . vs. 0^ for a step input of .80. 
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3 . Physical analysis of Third Order System with Backlash and Load 
Inertia* 

The backlash existing in the gear train between motor and load is 
placed inside the feedback loop as illustrated in Fig* 12. I'Jhen the 
backlash is taken up, the differential equations for the system are: 



i/n ( * i/n is) 



T 



C 




( 20 ) 




I + PI 



a 



a 



( 21 ) 



E 



- h<>c 



( 22 ) 



Substitute equation (22) into equation (21) it becomes: 




‘‘aVc - I, * 



a 



(23) 



Differentiate equation (19) we obtain: 



Tc » t + 1/H (K^f^, + 



(24) 



and 



I 



a 



l/K^(T^) 



(25) 





( 26 ) 
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Rearranging 


equation (26) 


and multiplied both sides by 


it becomes t 




2 


•• • 


2 


2 


K 1 K 2 K 30 R = 


1/N(N J + 
M 




1/N(N + f, + 

M L 


'"’■’m ^ ''V®c 




P/N(N^f„ + 
M 






(27) 
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Fig, 12. A block diagram of a third order servo with backlash 
and load inertia. 
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Where : 



= amplifier constant Volts/Volts 

ss motor torque constant ft-lbs/amp. 

K 2 = error measurement constant Volts/rad, 

K. = K /L = Volts/Volts-henry 

J a 

P = R/L = Armature resistance = Ohras/henry 

Armature inductance 

R = Armature resistance = Ohms 

L = Armature inductance == henry 

N = Gear ration ^ = rad. motor/rad. output 

2 

J.. = Inertia of system without load ft-lbs-sec /rad 

M 

2 

J * Inertia of load ft-lbs-sec /rad 

Li 

f . = Friction of system without load 

M 

f Friction of load only 

J = + J, K = K„K„ 

ML 16 

•G = Displacement of the combined system. 

6 = Displacement of the load only. 

Li 

6 s Input displacement. 

K 
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Letting N » Ij equation (27) becomes: 



( 28 ) 



This is a third order differential equation, hence we can not solve 
it by using phase plane technique. In order to be able to draw the phase 
plane trajectory we have to solve the above differential equation either 
by longhand or by using analog computer. Solving equation (28) as a fun- 
ction of time (t). 






yjK, K 



9 . 



R + 






^ c-f/j + Bjs\(p-r/j)s 

C f/J t + C P-f/ + ( K, K/j)J 

^ Cf/J -f 

CP/J+PJS+ c-f-P/j-Js -f ck,k/j-j 



letting 

- 
C ~ 



K^K 



1 and 0. 



1/s, equation (29) becomes: 



1/JK.K t[^-r 0/JtP;.S% (Pf/J)SJO ,, 
^ S-f (-f/J+ P)J ■»* 



LO 



= 

or c 



+ 



c -f/j t PJ -f (T pp/sjs 4 ‘/j 

cvJt?)s''+ 

S[("s+aj(fs + b+jcJ(T- 5 >' -i-b -jc)] 

rs + C^/J--rP;34o ^co 

(s -r a J c s + b + j c J ^ -r b -JC ) 



(29) 



(30) 



(31) 



where i 



f/j + P 

4 P/ J 



c\ -p zh 
Zc\b ^ bV*" 
c\ C + c 



14 it^ilial CO if 2 cj it I o 12 ,s = O ^ ec^u io \2 i 50 becomes 



0 - ^ 
" J 



_ CA 



_(b~jc)t 



o C b^ ‘t~ c ) 



+ 



^ ZjC(a-b-jc)' 
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+ e 



(b - jc)t 



2jc(a - b + jc) 



u(t) 



(32) 






i/j 



(33) 



+ a) (s + b + jc) (s + b 



jc) 



e. 



-at 



-( b - jc)t 



•f e 



(b - jc)t 



2 2 2 
-2ab + b c 



u(t) 



2jc (a - b - jc) 2jc( a - b + jc) 



(34) 



From equations (32) and (34) we can draw the phase plane trajectory 
of the velocity vSp displacement (see Figs 13) 







(1/J)s 

(s-fa) (s+b-^ jc) ( s + b - jc) 



(35) 
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33 



^ I 

i 




^-(b - jc)c 
2jc(a - b - jc.) 



u(t) 



2jc{a - b> - Jc) 



(36) 



Fig. l'^ ifj the plot of scceleration vs. displaic6f!.ient of th* combined 
system. 

The equ*'ntion fot. the drifting load Is detetriiined by the friction and 
inertia of the himd. 






(37) 






(38) 



L L 

This slope detertfiin/>s the point of seperation. 

Solving eqnat ion (i7' as a function of time (t)g, it becoTrfis; 



I 0, + 1 e, 

Lo — y y — i~ Lo 

■-{ s(s + d ) 



(39) 



wheic 



u 0^, - imitial comdition of displacement 

l.D Co Mo 



9^ 5is initial ccndltiori of velocity 

Lo (:o Mo 
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IS 




(40, 



[ 





A 



-dr 




i!!(t } 



o 





I 

( s + d ) 




{4!) 








(42) 



After the load has beeB sep«rered_, the equatiorj of the motor becomes; 









(43) 



J_ <5,, + l_ 

s s f s + d) 



Mo 



(44) 



Substitute this equation (44) loto equation (43) and solve as a function 
of a s and lettlr.g K.K - j' r becomes; 



1/J f 1 - 0 ) 

' M ' Me/ 



^ (1 + P)s + kPT 



+ ®Mc 



M Me 

^ S { £ 4 d) (s^ -r- { V f P)s + kP) 



s + k + P) { s -f d) &, 



M.O 



'Mo 



( S 4* ( k 4. p) S + kP) 



(45) 



where 






36 



nr 









Jj^s ( s + k) ( s + P) 



Mo 






^\~ 1 + J <s + k + ?)(£ + d 
+ rtJ- ” 

J s(s “ d) (s - k> {s - P) 



0, 



Mo 



M 



s (s f k) ( 8 + P) 



^ ^ ~ 

Jj^s(s + k) < S 4 P) 



e 



Mo + 



(46) 

-1 + ( s 4 b 4 P)(s 4 



Mo 



h 



e 

4 Mo 



(s 4 k)(s + P) 



Jj^( S = d) ( S 4 k)( s 4 P) 



(47) 



9, 



M = 



0 ■ «Ko> 



4 Mo 



[- 7 4 J 

4 L M 



( 8 4 k 4 P) (s 



4 d)]?' 



Mo 



J l(s 4 



4 P) 



J, 



M 



Jj_^(s 4 d) (s 4 k) (s 4 P) 



S0. 



Mo 



(48) 



^ k) ( S 4 P) 



'M 



_1__ i _ ^ 

^M t Ik 



- ( k 4 P) 4 P e - k e 

(Ik)^ (Pk)^( P - k) 



-1 



®Mo 



1 - t 



(k 4 P) 4 - P^e 



M 



+ Mo 



Pk (Pk)' 

ni<i<k + p) . 



M 



dkP 



®Mo 


4 


•'m 





fi * j^(k - aj] 



-kt 



k(P - k) ( k - d) 



(Pk) ( P-k ) 

-dt A 



d(k - d) (P - d) 



1 Jm<p - <<) - 1 ] I. 

P(P - d) (p - k) J 
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e 



M 









+ "®Mo 



M 



1 


4 |,ke'^‘^-Pe"*^*’’ . 


L 11 


u(t) ( 49 ) 


\Pk 


Pk (P - k) 


j. 




1 


-kt 

- e + 


-Pt ] 

e t 4 


“kP 


k(P - k) 


P(P - 


1^) J 


1 


1 


-Pt 

e 


u 


kP 


k(P - k) 


P(P - 


Mj 



Mo 



, -kt , -Ft. 

-ke ; 



^ [■ (P - k)e + (P - d)e ^*"+ (k - d)e 
(k - d) (P - d) (P - k) 



e 



Mo 



M 



_ ^ -kt "Pt. 
-g ) 

(P - k) 



u(t) 



( 50 ) 



e 



M 



22 


1 

f 


r -kc -Pt 
1 e -4 


1 

oi 

® i 

1 

1 

1 








\ )P - k) 


j 


\(P - k) J 



K (-^mC +Pe’^^) 

+ Mo ) M' -^ + 



'M 



(P - k) 



[d(P - k)e'*^'‘^-k(P-d)e'*^’^+P(k - d)e'^*^]|^ 
(k - d) (P - k) (P - d) f 



rt f , -kt „ -Pt 
4 0 .. j -ke 4 Pe 
Mo i 



(P - k) 



u(t) 



( 51 ) 
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Flgo 15 is the plot of velocity vSo displacemeiat of the motor and 
the plot cf displacement of load vSo displacement of the motor in the 
backlash regiono 

From the curve of the load-displacement vs* the motor-displacement 
we can determine the point of recombination for the motor* Knowing 
the displacement of the motorj the velocity of the motor at recombina- 
tion can also be determined* 

Applied the law of conservation of mamentums 






M 






Vm 



(52) 



From the above equation 0 and 0 are known quantities, so we 

JL n 

can evaluate the value of 0^o 

For the determination of the initial condition of the velocity at 
the point of recombination we have to apply the torque-equation of the 
motor* 



T 



M 



J. 



MM 




(53) 



Multiply equation (51) by and equation (50) by and add 

it together we obtain? 



T 



M 




-Pt 







Mo 



-kt -Pt 
\ t e -e 

I (P - k) 
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+ 6 ,, 

Mo M 



-kt -pt -dt 

(ke +Pe > + d(P - k)e -k(P 



(P - k) 



d)e'^^+ P(k - + 

J^(k - d) (P - d) (P - k) J 









1 - ke'^Vpe~^^ 

I (P - k) 




, „ -kt „ -Pt, ) 
(-Pe +Ke ) i + 

(P - k)P J 



e. 



Mo 



/D "kt , -Pt. 
(Pe -ke ) 

(P - k)P 






k(Pe~^*^-ke~^*^) + C- k(P - k)e~^'^+k(P - d)e'^*^+k(k - d)e~^ ^ + 



(P - k) 



J^(k - d) (P 



d) (P -k) 



J 



e J 

Mo M 






(P - k) 



-] 



u(t) 



(54) 



Figp 16 is the plot of the motor- torque vs. the motor-displacement. 
Knowing the motor-torque and the motor velocity at the point of re- 
comb inationj, we can determine the motor-acceleration at that point from 
equation (33), 

From here on we wolve for equation (32) using the initial condition 

3 * until it hit the point of separation again, 

Co Co Co 
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\ 
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43 
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Case Is 



For the coTblned syiten 



C » 1/ 

a 

c 


. i.s«^ * ic ^ °c ’■ 


1.0 








(55> 


Wheie 


J • 1.0 


£ = 


1.0 


N =: 


1.0 






J„ - -6 




.68 


P = 


1.0 






Jl .. .4 




,32 


K^K « 


1.0 














1.0 




G ( s ) » 


i 






1 






s ( + 1.8s 4 1.0) 




s(s + 


.9 + j.436) 


(s -f 


,9 - .J436) 



Fig. 17 shews the root locus of the system to be irsvestigaced. To 
solve equation (557 , useo an analog computer and the result is plotted 
in Fig? (18). 
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Th-2 load fiepa»:at(?d with the 'Squation? 



,40. 



, 326 ^ 0 

L 



The slope of the load is 






-.80 



(57) 



Eqiii»tloTi« (56) is plotted In Fig, (18). 
The fqvjatton of the niccor after separationr 



+ 7,467*>jj t <■ 1,6670^^ « 1.667 (58) 



Fro!f« tfts slope of the load v/a can find the initial condition for 
equation (58), In our case these are; 





« e 

l-o 


“ e,. 

Ho — 


1 . 360 


®c. 


1<jO 


a 0 

Mo 


.424 


v. 


3 ■ 

SC 0- 

Lo 


^Mo * 


,324 



lislng fhesf inttiai conditions we solve for equar.icn (58) and 

the reS'oJt is plotted in Fig, (IS), 

From Fig, (23) for a backlash » .42 rad. 

0. = -.074 rad/ sec, 

L 

0., “ .748 rad 0: ™ .480 rad/sec. 



8^ .328 rad 

L 



M 



M 



.748 rad 
,450 ft -lbs 
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o6(.480) 



.4(,074) 



,258 rad/ sec. 



0^ -- - .258 



192 rad/sec , 



From Fig. (28), for 0^ » .258 rad/sec, 0^ = .328 rad. 

Initial conditions are; 

0„ » . 328 rad. 

Co 

6„ s ,258 rad/ sec, 

Co 

0^^ = ,192 rad/sec^. 

Using these as Initial conditions we solve for equation (55) 

* *• * 

agaioc Values of 6 _ 0 and 0* at points of separation are tabu- 

c c c 

lated Dll table (5), with several values of backlash. 
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TABLE 5 







% 


Backlash 


1.360 


.424 


-.324 


1,0 


.600 


-.430 


.336 


1.0 


1.377 


.444 


-.342 


1.0 


1,360 


.424 


-.324 


.9 


.620 


-.420 


.326 


.9 


1.359 


,420 


-.320 


.9 


1.360 


,424 


-.324 


.8 


,639 


-.428 


.308 


.8 


1.348 


.408 


-.320 


.8 


. 666 


-.396 


.294 


.8 


1.321 


.380 


-.300 


.8 


o 

o 


-.390 


o 

CM 

O 


.8 


1.320 


.376 


-.298 


.8 


1.360 


.424 


-.324 


.7 


.652 


CM 

0 

1 


.310 


.7 


1.332 


,392 


-.304 


.7 



I 
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QD 



TABLE 5 continued 







» 


Backlash 


1.360 


.424 


-.324 


.6 


.660 


-.396 


.292 


.6 


1.306 


.356 


-.286 


.6 


.728 


-.316 


.224 


.6 


1.256 


.300 


-.234 


.6 


1.360 


.424 


-.324 


.5 


.666 


-.394 


.288 


.5 


1.286 


.334 


-.270 


.5 


.736 


-.312 


.226 


.5 


1.240 


.280 


-.222 


.5 


.768 


-.276 


.192 


.5 


1.220 


,254 


-.204 


.5 


.788 


-.250 


.186 


.5 


1.200 


.230 


-.200 


.5 


.800 


-.240 


.170 


.5 


1.192 


.220 


-.188 


.5 


.810 


-.230 


.165 


.5 
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TABLE 5 continued 





o 


'^c 


Backlash 


1.360 


.424 


-.324 


.4 


.676 


-.382 


.264 


.4 


1.260 


.306 


-.244 


.4 


.754 


-.290 


.200 


.4 


1.208 


.246 


-.194 


.4 


.798 


- . 240 


.168 


.4 


1.164 


.216 


-.166 


.4 


.825 


-.208 


.144 


.4 


1.140 


.200 


-.142 


.4 


.832 


-.196 


.118 


.4 


1.120 


.194 


-.134 


.4 


1.360 


,424 


-.324 


.3 


.705 


-.344 


.264 


.3 


1.240 


.280 


-.234 


.3 


.800 


-.234 


.172 


.3 


1.175 


.208 


-.160 


.3 


.846 


-.180 


.134 


.3 


1.155 


.184 


-.140 


.3 


.852 


-.170 


.120 


.3 


1.132 


.158 


-.120 


.3 


.864 


-.158 


.110 


.3 
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rx\ 



I 



TABLE 5 conttnued 











Backlash 


1.360 


.424 




-.324 


.2 


.720 


-.344 




.254 


.2 


1.222 


.276 




-.210 


.2 


.820 


-.220 




.198 


.2 


1.136 


,170 




-.120 


.2 


.880 


-.148 




.108 


.2 


1.098 


.124 




-.100 


.2 


1.360 


,424 




1 

1 

1 

1 1 

L. 1 

1 

1 

f 


.1 


.740 


-.312 




,238 


.1 


.1200 


.236 




1 

• 

CD 

O 


.1 


.860 


-.190 




.140 


.1 


1.120 


.134 




-.112 


.1 


.910 


-.110 




.094 


.1 


1.074 


.092 




-.064 


.1 


.940 


-,080 




.052 


.1 
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Values of 6., 0, , 0..„ Q.,. 0_, 0_ and T,. at points of recosnbina- 

L L M M O L M 

tion are tabulated on Table 6 for several values of backlash. 



TABLE 6 











T 

M 






Backlash 


1.804 


.070 


.804 


-.842 


-.615 


-.477 


-.233 


1.0 


.124 


-.050 


1.124 


.904 


,632 


.522 


.214 


1,0 


1.840 


.076 


.840 


-.852 


-.636 


-.481 


-.251 


1.0 


1.798 


.074 


.898 


-.800 


-.604 


-.4504 


-.244 


.9 


.170 


-.060 


1.700 


.790 


.606 


.450 


.246 


,9 


1.794 


.670 


.894 


-.788 


-.590 


-.445 


-.234 


.9 


1.790 


,080 


.990 


-.750 


-.596 


-.418 


-.262 


.8 


.86 


-.072 


.986 


.746 


,584 


.423 


.246 


.8 


1.784 


.082 


1.084 


-.692 


-.582 


-.383 


-.275 


.7 


.206 


-.056 


.906 


,664 


.560 


,376 


259 


.7 


1.734 


.072 


1.034 


- . 680 


-.548 


-,379 


-.245 


.7 


1.776 


,090 


1.176 


-.632 


-.550 


-.343 


-.280 


.6 


.258 


-.072 


,858 


.610 


.526 


.3372 


.2562 


.6 


1.670 


.068 


1.070 


-.376 


-.350 


-.1984 


-.1912 


.6 


.398 


-.054 


,998 


.554 


.440 


.3108 


.1914 


.6 


1.558 


.058 


.958 


-.539 


-.422 


-.2996 


-.1824 


.6 


1.762 


.104 


1.262 


-.560 


-.552 


-.294 


-.316 


.5 


.274 


-.080 


.774 


.544 


.494 


.2944 


.2584 


.5 


1.602 


.082 


1.102 


-.512 


-.448 


-.2744 


-.2286 


.5 


.404 


-.046 


.904 


.488 


.406 


.2744 


.1864 


.5 


1.508 


.066 


1.008 


-.484 


-.390 


-.264 


-.1788 


.5 


.472 


-.042 


.972 


.460 


,370 


.2592 


.1626 


.5 


1.472 


.055 


972 


-.476 


-.364 


-.2636 


-.1332 


.5 


.512 


-.034 


1.012 


.448 


.334 


.2552 


.1298 


.5 


1.430 


.048 


.930 


-.434 


-.324 


-.2412 


-.131 


.5 


.536 


-.032 


1.036 


.422 


.328 


.2404 


.1376 


.5 


1.414 


.041 


.914 


-.426 


-.310 


-.2364 


-.1208 


.5 


Figs. 


18, 19, 


20, 21, 


22, 23, 24, 25, 26 


j, and 27 show us the 


result of 


the combined 


trajectories with 


backlash 


from ol till l«0cie^,ln our case 


limit 


cycle occurSo 


Values of 


limit cycles versus 


backlash 


is tabulated 



on Table 7 . From the table we come to the conclusion that, if we in- 
crease the amount of backlash the amplitude of the limit cycle will in 
crease too. 
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TABLE 7 



For Gain-constant is l.O 



Backlash 


Amplitude of limit cycle 


.1 

.2 

.3 


.355 


.5 


.435 


.6 


.5625 


.7 


.725 


.8 


.760 


.9 


.810 


1.0 


.840 
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^0 Aiaalcg ccmpvrcer <» 

Foe th.ef combined sy'£'rems 



\ 



(59) 



This i« An equation of a siUM).er with 4 inputs loec 

? f «i> 

O^s, 6^ ani 






- 0 , 



0, 






Oc] 



fl 



(60) 

(61) 

(62) 



Equations (fcO), (61) and (62) are equations for an integrator. 
Hence for the combined system we need 1 summer 3 irtegrators and in 



addition 3 sign charge* s. 



Summary for the combined systemo 

The output of amplifier 1 which in our case is & is fed for- 

ward into amplifier 2 ( an integrator). One output of this amplifier 

is fed back into amplifier 1 through a sign changer amplifier 6. 

The other output is fed forward into another integrator i.e. amplifier 3 

Output #i of amplifier 3 is fed back into amplifier 1, 

Output #2 is fed forward Into integrator amplifier 4. The output of this 

is fed back through a sign changer into amplifier 1. The output of 

amplifier 2 and amplifier 4 are recorded into a X - Y plotter i.e. 

the output of amplifier 2 i.e. 6 Is fed into the Y - axis and the out- 

u 

put of amplifier 4 I.e. 8 is fed into the X - axis of the X-Y plotter. 

u 

So we have a recording of acceleration versus displacement of the com- 
bined systenic 

The other recording is output of amplifier 7, which is 0 is 

u 

fed into the Y - axis and the output of amplifier 4j which is 0 
is fed into the X - axis of the X-Y plotter. From this recording 
v;e have the phase - plane trajectory of velocity versus displacement of 
the combined system 

The circuit diagram of the analog computer set-up of the combined 
system is shown in Flgo 30. 



66 



Ftg. 30, A circuit diagram of computer set up # 1 for the combined system. 
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Calculation for the computer set up # 1. 

Amplifier l o 

W, s ^ 1.0 a, s 1.0 



W. 



R-. 



1.0 



= 1®0 



w. 



®3‘*^fl 



1.0 



a^ ” l.Ox.3 ~ .300 



W. 



V'fi 



1.8 



» 1.8 X .3 = .540 



Amplifier 2 . 



W, 



m 





1.0 



Amplifier 3. 



W. 




R. 



1.0 



A.mplifier 4. 



W, 



<=fV'7 



1.0 



fig. » loO 



a. * loO 
6 



Bj » loO 
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Amplifier 5 



Amplifier 6, 



Amplifier 7 . 



W. 



®8°^f5 



1.0 










1.0 



^lO ^ 7 = 1.0 

^0 



s 1,0 



Sig » 1.0 



a 



10 



1.0 



The components used in this*' computer set up # 1 are: 



"l 


S 


1.0 


M-Ohms 


Rfl 




o 

o 

?-4 


Meg. 


^2 




1.0 


Meg 










"3 


- 


.3 


Meg 












2 


.4 


Meg 












2 


1.0 


Meg 


CM 


= 


1.0 


microfarad 






1.0 


Meg 


"f3 




1.0 


microfarad 


h 




1.0 


Meg 


"f4 


« 
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The time scaling loO sec. 

The magnitude scaling for the displacements 50 Volts = 1,0 rad. 

The magnitude scaling for the velocity s 50 Volts * 1,0 rad/sec, 

2 

The magnitude scaling for the accelerations 50 Volts =1,0 rad/sec , 



When the system separate,, the equation of the load becomes s 




\ M 

■ -i- [ 

The equation of the motor after separation becomes s 




(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 
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Equation (64) and (65) are the equatioi^s for Integrators, 

For the load equation we need two Integrators and one sign changer. 

For the motor we need 3 Integrators, one summer and 2 sign changers. 

The summery of the computer set up if 2 Is as follows: 

One of the output of amplifier 1, which Is 0^ is fed back Into amplifier 
1. The other one Is fed forward into amplifier 2, This amplifier has an 
output of (load displacement). The output of amplifier 3 is fed 

into tl$a Y-axis of a X-Y plotter, while the X~ axis gats Its Input 

from amplifier 7. So we have a recording of load displacement vs. 

motor displacement. 

Amplifier 4 has 4 Inputs l,e. 0^^, 0^^, 0^^ and 0^^. One of the 

output of amplifier 4 is fed back into amplifier 4. The other one is 

fed forward Into amplifier 3. There are two outputs of amplifier 5, 
one Is fed back Into amplifier 4 through a sign changer amplifier 8. 

The other one Is fed forward Into amplifier 6. Amplifier 9 has 2 Inputs 
t.e. Oj^ and 0^^. 

Recording points are: 

0 from amplifier 5 is fed Into the Y “ axis of a X-Y plotter. 

0.. from amplifier 7 is fed into the X - axis of the X-Y plotter. 
M 

Hence we have a phase plane trajectory of motor ■velocity versus motor dis- 
placement. 

Motor torque T^^ from amplifier 9 is fed into the Y - axis and 
the motor displacement from amplifier 7 is fed Into the X - 
axis of the X-Y plotter. 

From this we have t!ie phase plane trajectory of motor torque versus motor 
displacement. The circuit diagram of the analog computer set up for the 
load and the motor after separation is shown in Fig. 31. 
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s circuit diagram of computer set up # 2 for the load and motor 
after separation* 
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Calculation for the computer set up # 2. 
Amplifier 1 . 



W 



1 



CfioR^ 



.800 



Amplifier 2, 



W. 



C ^ 2*^2 



1.0 



Amplifier 3 . 



W, 



®3'^f3 



1.0 



Amplifier 4 . 



Amplifier 6 , 



W, 



W. 



W, 



W. 



Amplifier 5. 



W, 












^f4’\ 



^f4'^7 



= .547 



.547 



.800 



s 1 o800 



Wq ^ 8 ® loO 

"f5-«8 



C^,.R„ 
f6 9 



1.0 



Amplifier 



W 



10 



®10'*^f7 



1.0 



R, 



10 



Amplifier 8 . 



W_ ® ^ 1.0 

11 " ' 



11 
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» .800 

- 1.0 

S 

- .547 

. .547 

a^ “ .800 

a^ ss 1.8 X .5005 » .901 

3g E 1.0 

= 1*0 



11 



» 1.0 



Amplifier 9 . 

® ^12° ^f9 s .480 a.„ » 

12 — 12 

^12 

Wj3 ® ^13*^f9 s .600 " 

*^13 



.480 



.600 
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The component ^sed in the compniter set up # 2 ares 
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The time scaling and magnitude scaling used in the computer set up 
# 2 are the same as used in the previous computer set up. 



5. The result of Investigation. 

We run experiments with different values of the system gain i.e. 

6) K « .95 

b) K * 1.00 

c) K * 1.20 

The result is plotted in Fig. 32, This figure shows the graph of 
the amplitude of limit cycle versus the backlash with different values 
of K. The graph tells us, that when we increase the amount of backlash, 
the amplitude of the limit cycle will increase too, provided we keep the 
gain constant. Now if we increase the gain, letting the amount of back- 
lash constant, the amplitude of the limit cycle will increase. Notice 
that by decreasing the amount of the gain by ,05 (l,e, from K => 1,0 
to K s .95), the curve moves more to the right as does it to the left 
if we increase the system gain from K a 1,0 to K s 1,20, The values 
of the limit cycle with different gain constant is tabulated on table 



8 . 



TABLE 8 



Backlash 

in 

rad. 


K - .95 


Amplitude of 
K s 1.00 


limit cycle in rad 
K s 1.20 


1.0 


.598 


.840 


-- 


.9 


.540 


.810 


— 


.8 


.425 


.760 


-- 


.7 


.305 


.725 


.970 


.6 


-- 


.5625 


.900 


.5 


-- 


.435 


.790 


.4 


-- 


.355 


.671 


.3 


-- 


-- 


.552 


.2 


-- 


-- 


— 


.1 


• • 


- - 


• — 
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3o Discussion and Conclusion 



It is seen from the results of this study that backlash in a third 
order system may cause a limit cycle or may not cause a limit cycle 
depending on the amount of backlash and the system gain. With large 
amounts of backlash and low gain a limit cycle is unlikely. As the 
gain is Increased it becomes more probable that a limit cycle will 
exist. 

If a limit cycle does exist the amplitude of the limit cycle depends 
on both the amount of backlash and the gain. At low gain the amplitude 
of the limit cycle (in radfans) may be considerably less than the 
amount of backlash (in radians). Conversely, at high gain the amplitude 
of the limit cycle may exceed the amount of backlash. 

These conditions have been determined from study of a system with 
adjustable gain, but without variation of the inertia and friction, nor 
of the amounts of inertia and friction on each side of the backlash. 



These factors need further investigation. 
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